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Abstract 

Bound vesicles subject to lateral forces such as arising from shear flow are 
investigated theoretically by combining a lubrication analysis of the bound 
part with a scaling approach to the global motion. A minor inclination of 
the bound part leads to significant lift due to the additive effects of lateral 
and tank-treading motions. With increasing shear rate, the vesicle unbinds 
from the substrate at a critical value. Estimates are in agreement with recent 
experimental data. 

PACS: 82.70-y, 87.45-k, 47.15 Gf. 

Introduction. The equilibrium aspects of the interactions between membranes or vesicles 
and substrates have been explored intensely over the last decade |J. Quantitative exper- 
imental data for both the mean shape and the fluctuations of the bound part of a vesicle 
have been obtained by using phase contrast microscopy and reflection interference contrast 
microscopy (RICM) [0. A qualitative new step concerns the study of bound membranes 
under controlled non- equilibrium conditions such as the behavior of bound vesicles under 
shear flow. Apart from its fundamental significance, this system can serve as a model for 
the biologically ubiquitous situation of adhesion of membranes under flow. A prominent 
example occurs for leucocyte or platelet adhesion in capillary flow. Clearly, for biological 
systems, the non-trivial kinetics of specific adhesion molecule pairs under a ramped force 
contributes essentially to the dynamic unbinding of these cells (see, e.g., [||-|5| and refer- 
ences therein). Still, a thorough understanding of the model case of a bound vesicle with 
its interplay between unspecific interactions and flow will be a prerequisite for gaining a 
comprehensive picture of these important dynamic interactions. 

A significant experimental step in this direction has been achieved recently by combining 
RICM with a flow chamber |5[]. With this set-up one can study the configurations of bound 
vesicles under shear flow. It was observed that these vesicles unbind from the substrate at 
a critical shear rate. However, the effective lift force was found to be about two orders of 
magnitude larger than what was predicted in previous theoretical work ||]. 

The purpose of this letter is to analyze theoretically the dynamically induced interaction 
between a substrate and a bound vesicle under a lateral force such as arising from shear flow. 
This problem is challenging since it involves two vastly different length scales. Typically, the 
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vesicle size is of the order of 10 /iin whereas the distance between substrate and vesicle is 
of order 10 nm. A brute force approach trying to solve numerically the equations of motion 
of such a configuration as it has been done for free vesicles in shear flow J7| is bound to 
require a very fine discretization and, consequently, to face high computational costs. A 
first step in this direction has been achieved recently for the computationally less expensive 
two-dimensional case [@.[|. 

For the experimentally relevant three-dimensional two step approach will be 

followed here. First, the bound part of the vesicle will be treated quantitatively within 
the lubrication approximation which holds if the lateral extension of the bound part is 
significantly larger than the distance from the substrate. As a result we will find that 
whenever this bound part is tilted a significant hydrodynamic lift arises even for small tilt 
due to the additive effects of translation and relative membrane motion, i.e. tank-treading. 
In a second step, we couple this lubrication analysis into a scaling approach of the overall 
vesicle motion. As a result, we predict a critical lateral force beyond which vesicles will 
detach from the bound state and, thus, undergo a dynamically induced unbinding transition. 
This transition must be distinguished from an equilibrium unbinding transition due either 



to fluctuations ||10| , pT| or a competition between adhesion energy and curvature energy |12 

Geometry. In equilibrium, a vesicle bound to a substrate by a potential V{h) acquires 
a spherical cap-like shape if the depth W of the potential is sufficiently deep, see Fig. 1. 
The shape can then be characterized by two parameters, the radius R of the spherical cap 
and the radius R a < R/2 of the adhesion disc. The distance ho of the adhesion disc from 
the substrate is determined by the location of the minimum of the adhesion potential. In 
equilibrium, the adhesion disc is parallel to the substrate but we will here allow a small 
inclination or tilt angle a. 

We now apply a force F x parallel to the substrate. The physical origin of this force 
can either be a linear shear field with shear rate 7 or a gradient in adhesion energy VW 
The force then scales as F x ~ 7?7-R 2 , or F x ~ VWR 2 ., respectively. As a result the 
vesicle moves with a velocity v into the same direction. Since the membrane is fluid, we have 
to allow for tank-treading motion which we assume for the spherical part to be a uniform 
rotation in the x, z plane at an angular speed v m /R. This membrane flow on the spherical 
part enters (or leaves) the rim of the adhesion disc with a y dependent velocity 

v m {y) = v m {l-y 2 /R 2 ) l l\ (1) 

where \y\ < R a . The velocities v and v m will later be determined from force balances but 
for the moment they are assumed to be given. We show first that such a motion generates 
a significant hydrodynamic lift on the vesicle. 

Lubrication theory. For future reference, the lubrication approximation will be set up not 
just for a tilted adhesion disc but for a general membrane configuration parametrized by h(r) 
in a Monge representation above a substrate at z = with r = (x, y). The lateral extension 
R a of the membrane is much larger than the average height ho above the substrate. The 
velocity field is written as v(r, z) = u(r, z)+w(r, z)e z where u(r, z) is the component parallel 
to the substrate. At the membrane, we specify the velocity as v(r, h(r)) = u h (r) + w h (r)e z . 
At the substrate, no-slip boundary conditions imply v(r, 0) = 0. 

The Stokes equations for the incompressible fluid between substrate and membrane read 
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V (d 2 + V 2 )u = Vp, (2) 



and 

r)(d 2 z + V 2 )w = d zP . (3) 

Here V is the gradient operator in the plane, i.e. V = (d x ,d y ) in Cartesian coordinates. 
From the continuity equation d z w + Vu = 0, it follows that w/\u\ ~ O(h /R a ) for small 
ho/R a . This scaling is the essential observation in lubrication theory. It implies via (0) 
and (H) that the pressure is a function of r but independent of z to leading order in h /R a . 
Likewise the V 2 term can be ignored compared to the d 2 term in (|2J). The latter equation 
can hence be integrated as 

u(r, z) = Vp(r) z{z - h)/2rj + u h {r)z/h, (4) 

which satisfies the boundary conditions at z = and z = h. Applying the V operator to this 
equation, inserting the resulting expression into the equation of continuity and integrating 
the latter over z from to h yields a Reynolds-type equation 

V 2 p + 3VpVh/h = 12r]w h /h 3 + 6rjV(u h /h)/h. (5) 

This equation for the pressure holds for any membrane configuration in the lubrication 
approximation. We now specialize to the circular adhesion disc of radius R a tilted at a small 
angle a, i.e. h(r) = h Q + ax, see Fig. 1. No slip boundary conditions between membrane and 
fluid imply for the translational and tank-treading motion introduced above the boundary 
values w h (r) = —av m (y) and u h (r) = {v — v m (y))e x for small a. With Vh = ae x , equation 
(|) becomes 

V 2 p(r) w -6ar](v m (y) + v)/h 3 (6) 

to lowest order in a. First, assume that there was no tank-treading motion, v m = 0, and, 
hence, no r dependence of the inhomogeneity in this Poisson equation. The solution then is 

p(r) =Zar ] v(R 2 a -r 2 )/2hl+p(R a )(r/R a ) C os<P (7) 

where p{R a ) cos <p is the ambient pressure along the rim of the adhesion disc parametrized by 
the azimuthal angle <fi. Integrating the excess pressure over the entire adhesion disc yields 
the total lift force on the membrane disc as 

Fi = J dA p(r) = 3Ttar]vR 4 Jih 3 . (8) 

For v m 7^ 0, the solution of the Poisson equation @ is slightly more involved because of the 
y dependence of the rhs. We are interested only in the total lift, which involves averaging 
over the whole adhesion disc. Given the form (II), it is then clear that (|8|) still holds with 
v replaced by v m up to a dimensionless function f(R a /R) of order unity which will not be 
needed for the scaling analysis to follow. The important point is that both translational 
motion and tank-treading motion contribute similarly to the hydrodynamic lift (|8]). This 
lift increases strongly with decreasing distance ho of the membrane from the substrate. Note 
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that the reversibility of the Stokes equations implies that there is no hydrodynamic lift for 
a non-tilted configuration with a = 0. 

Scaling analysis of global motion. We now have to link this lubrication analysis of the 
adhesion disc to the overall vesicle dynamics. We apply the force F x parallel to the substrate 
and assume first the rotational degree of freedom being locked at a = 0. Then, there 
are two conditions which fix the two velocity parameters v and v m uniquely. We write 
these conditions in scaling form which means that we ignore all numerical prefactors of 
the respective terms with the understanding that these prefactors do not exhibit singular 
behavior for small h /R a . 

(i) Force balance in the x direction: 

F^rjAvRl/ho + rjvR, (9) 

where Av = v — v m is the velocity of the bound part of the vesicle relative to the substrate. 
The first term on the rhs is the lateral force exerted by the lubrication layer, the second 
term is the hydrodynamic drag of the exterior fluid outside of the lubrication layer. 

(ii) Dissipation balance: 

F x v ~ r](Av) 2 Rl/h + V v 2 R + nv 2 m R. (10) 

The lhs represents the power applied by the external force on the system. The first term 
on the rhs is the dissipation in the lubrication layer. The second one is dissipation in the 
exterior fluid outside of the lubrication layer. The third term is the dissipation within the 
vesicle due to tank-treading. We now divide the first equation by rjR, the second one by 
rjvR and introduce with vs = F x /6ttt]R the Stokes velocity of a spherical vesicle in infinite 
space. 

The two equations (^],|TUJ) then read 



v s ~(3Av + v (11) 

and 

v s ~ (3{Av) 2 /v + v + v 2 Jv (12) 

where (3 = R 2 /h Q R is a dimensionless variable measuring the relevance of the substrate. For 
the scaling analysis, we have to distinguish three cases. 

(i) For (3 « 1, the vesicle is too far away from the substrate to be affected significantly 
in its motion. We will not consider this case further. 

(ii) For (3 ~ 1, we have Av ~ v ~ v m ~ v$- In this case, both translational and 
tank-treading velocity are of the order of the Stokes velocity. 

(iii) For f3 » 1, equations ([TT|JT2]) imply Av ~ 1/(3 and, consequently, v m ~ v ~ v$- In 



this limit, tank-treading and translational velocity become equal and both are of the order 
of the Stokes velocity 0. Tanktreading thus restores (up to factors of order unity) the free 
mobility which would be impossible for a rigid object with finite R a so close to a substrate. 
In summary, we can write in each of the two interesting cases 

v s ~ F x /r}R and Av/v ~ h R/R 2 a . (13) 
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Hydrodynamic lift. Using these results, which will hold within perturbation theory also 
for small nonzero tilt angle a, the lift @ can be written as 



F l ~ a V vR A Jhl ~ aF x (R a /R)(R a /h ) 3 . (14) 

This expression still depends on the unknown tilt angle a which must be determined next. 

In general, the transversal motion considered so far for the rotationally locked shape at 
a = generates a torque M acting in the x, z plane. The origin of this torque are the 
lateral force and hydrodynamic interactions. Their sum will scale as M ~ rjvR 2 ~ F X R. 
Counteracting to such a torque is a torque arising from the confining adhesion potential 
which favors a = 0. The energy E(a) of a tilted adhesion disc compared to a non-tilted one is 
given by E(a) ~ (ir/4)a 2 V"R^ where V" = V"(h ) is the curvature of the adhesion potential 
at the minimum ho- Balancing the torque d a E derived herefrom with the hydrodynamic 
one leads to r/vR 2 ~ aV"R\ or 

a ~ nvR 2 /V"R 4 a ~ F x R/V"R 4 a . (15) 

We have written these relations as if a was positive. Thus, we have implicitly assumed that 
the hydrodynamic torque acts to increase a. While this cannot be proven without a full 
hydrodynamic calculation, evidence for this assumption in the case of shear flow arises from 
the fact that for free vesicles the strain component of the shear orients slightly elongated 
shapes towards a positive tilt angle f7j. Note that if the sign of a was reversed, the scaling 
relations derived so far would still hold with lift being replaced by an additional force pushing 
the vesicle towards the substrate. 

With this expression for the tilt angle a > 0, we can now calculate the lift force from 

as 



F t ~ tj v R / h^V" ~ F 2 /h 3 Q V". (16) 

For further evaluation, we need the specific form of the adhesion potential V(h). A fairly 
universal relation can be derived if we assume that this potential can be characterized by 
two scales only, the location h of its minimum and its depth W = \V(ho)\, see Fig. 1. 
Then, one has V" ~ W/h^, which implies F t ~ F 2 /Wh . 

For a comparison with experimental quantities, it is convenient to express the adhesion 
energy W in terms of an effective tension S using a Young Dupre equation |12||. Except in 



the tense spherical limit R a — > 0, both quantities are of the same order, i.e. £ ~ W . For 
vesicles in shear flow, we can thus write for the lift 

Fi ~ F 2 /Wh ~ (r / 2 7 2 J R 3 /S)(i?// i0 ). (17) 

This result is a factor of order R/ho ~ 10 2 — 10 3 larger than a previous theoretical estimate 
IH. On the basis of the latter, it was argued in Ref. 0] that the experimentally observed lift 
was 2 to 3 orders of magnitude larger than theoretically expected. In the light of the present 
theory, this apparent discrepancy is most likely due to the factor R/h missed in Ref. |4j]. 

Dynamical unbinding. The lift, if small enough, will displace the vesicle slightly from 
the static equilibrium at ho. The new dynamical equilibrium position can now be found by 
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balancing the lift with the restoring force, F z rs —TtR 2 a {h — ho)V", arising from an expansion 
of the potential around its minimum. The relative shift thus becomes 

(h - h )/h ~ F 2 JhtV" 2 Rl ~ F 2 JW 2 RI (18) 

Typically, if (h — h Q )/h ) ~ 1, the lift will be too strong to be compensated by the 
attractive potential. Then, the vesicle will unbind from the substrate under the action 
of a lateral force. Using this criterion, we can determine the critical lateral force F£ as 
F£ ~ h\y"R a ~ WR a . If the lateral force arises from a shear field, the critical shear rate is 

7 C ~ WR a /riR 2 . (19) 

Depending on the specific conditions, measurements for the adhesion energy W have 
obtained a vast range of values from 10 _1 to 10 -6 erg/cm 2 [pf , P , Z| . Consequently, the 



critical shear rate will also depend strongly on the conditions. Rather than using (|19|) to 
estimate this force, we chose a typical shear rate of 7 = 1/s and determine the critical 
adhesion energy W c at which the dynamical unbinding should occur. For a vesicle with 
R = 10/iin, R a /R = 0.1 and a fluid viscosity of rj = 10 _2 erg s/cm 3 , we find a critical 
adhesion energy of W c ~ 10~ 4 erg/cm 2 which is well within the above range. For W > W c , 
vesicles with the above geometry will unbind from the substrate. 
The tilt angle a c at the critical value is 

a c ~ Fl R/V"Ri ~ Fl Rh 2 /WR 4 a ~ h\ R/R 3 a . (20) 

Within this simple assumption about the adhesion potential, the scaling of the critical tilt 
is thus predicted to be determined exclusively by geometrical quantities. Note that since 
ho << R a <, R, the critical tilt angle can be very small, e. g., a c ~ 10 -3 for h = 10 nm 
and R, R a as above. 

Beyond unbinding. What happens with the vesicle for larger lateral forces F x > F£ ? If 
the adhesion potential decays to for large h, one can expect that the vesicle will continously 
drift away from the substrate presumably in a combination of tank-treading and tumbling 



motion. This regime has been studied perturbatively for large h /R in Ref. |T5 |. 

If, however, as often in experiments, the vesicle is filled with a slightly denser fluid than 
the surroundings, gravity will keep it close to the substrate since the lift necessarily gets 
weaker with increasing distance. Such a state where after detachment the vesicle still keeps 
its spherical cap like configuration and translates at about 100 nm above the substrate has 
been reported experimentally || . We close with a somewhat speculative theoretical analysis 
of this state. The jump in separation from the substrate of about 100 nm could indicate that 
a short range adhesion potential is no longer responsible to the energetics after detachment. 
Neither can the potential then exert a restoring torque. Therefore, the total hydrodynamic 
torque must vanish. At which angle a this happens (if at all) can only be determined by 
a full hydrodynamic theory. The present scaling approach, however, allows to relate this 
tilt angle a to the separation h by balancing the lift force ( ]T4| ) with the gravitational force 
F g = gApV, where g ~ 1000 cm/s 2 and Ap is the density difference. Such a balance 
leads to a ~ (h / R a ) 3 (F g /F x )(R/ R a ). Using typical experimental values F g ~ 10~ 7 erg/cm, 
F x ~ 10~ 8 erg/cm, h /R a — 0.1, and R/R a ~ 10, one finds with a ~ 0.1 that a quite 
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small tilt angle could generate enough lift to sustain such a stationary state |L6]|. It must 
be emphasized that this estimate is speculative in two respects. First, we have not shown 
that the full hydrodynamics allows such a torqueless state but have rather postulated it. 
Secondly, numerical prefactors could easily conspire and change such an estimate about one 
order of magnitude. 

Summarizing perspective. We have analysed the hydrodynamics of a bound vesicle under 
a lateral force allowing for tank-treading motion. A positive tilt of the bound part will 
generate lift. By combining a lubrication analysis with scaling arguments, a dynamically 
induced unbinding transition is predicted when the lateral force exceeds a critical value. Even 
though our analysis used for simplicity the spherical cap configuration, the scaling relations 
will also hold for weakly bound states with rounded contact regions ||12|| . More quantitative 
theoretical calculations will be needed to fill in numerical prefactors necessarily missing in 
such a scaling analysis. Such full scale calculations will also reveal the role of transients 
and assess whether the theoretically somewhat speculative stationary translating state after 
unbinding exists. On the experimental side, special efforts should now be undertaken to 
measure the here theoretically predicted but so far not yet experimentally reported tilt 
angle which, even if minute, is ultimately the source of the lift. 
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FIG. 1. A bound vesicle under a lateral force F x (left) in a potential well V(h) with minimum 
at ho and depth W (right). The radius of the adhesion disc is R a , its tilt angle is a. The vesicle 
translates at velocity v and tank-treads at velocity v m . This motion generates a lift F\. 
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